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Abstract 

The gauge-gravity duality relates on-shell graviton multi-point func- 
tions in an anti-de Sitter black brane background to stress-tensor cor- 
relation functions of strongly coupled fluids. We show that, for con- 
formal fluids whose gravitational dual is Einstein gravity, such cor- 
relation functions are universal in a particular kinematic region and 
depend only on the Mandelstam variable s. Thus, we extend the uni- 
versality of the two-point function or, equivalently, the universality of 
the shear viscosity to entropy ratio to multi-point correlation func- 
tions. We also consider the leading corrections to Einstein gravity in 
string theory. The effect of these corrections is similarly universal, 
inducing a specific difference in the angular dependence of the stress- 
tensor correlation functions in the same kinematic region. We discuss 
schematically how multi-particle correlations in heavy-ion collisions 
can be used to test these findings and, thereby, probe the existence 
of a gravitational dual to the quark-gluon plasma and its universality 
class. We also propose a new way to measure the shear viscosity to 
entropy ratio using multi-point correlation functions. 
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1 Introduction 



The gauge-gravity duality states that strongly coupled gauge theories in 
four dimensions have a dual description in terms of weakly coupled gravity 
theories in five- dimensional anti-de Sitter (5D AdS) space [TJ [2]. In some 
sense, the gauge theory lives at the outer boundary of the AdS spacetime, 
and so this is a holographic correspondence [3]. 

If the bulk gravity theory is Einstein's, then all of its on-shell amplitudes 
depend on a single dimensionful parameter, the 5D Newton's constant G§. 
Hence, appropriately chosen ratios of amplitudes do not depend on G$. For 
example, the ratio of the shear viscosity to the entropy density, rj/s = 1/4tt 
11], is related to a ratio of two-point functions of gravitons [5] and thus 
independent of G5. In string theory, the leading corrections to Einstein 
gravity are universal and depend on one additional dimensionful parameter, 
the string tension a'. 

So far, the emphasis in theoretical and experimental studies has been on 
the two-point functions; the shear viscosity, conductivity, etcetera. From the 
experimental side, the prototypical example of a strongly coupled fluid is the 
quark-gluon plasma (QGP) [6J. The QGP is produced in heavy-ion collisions 
in the Large Hadron Collider (LHC) at CERN and the Relativistic Heavy 
Ion Collider (RHIC) at Brookhaven. Experiments on the QGP have already 
had significant success at testing the duality. In particular, the unusually low 
value of r]/s was first predicted via the gauge-gravity duality [1] and then 
later substantiated [TJIH]. 

As multi-particle correlations in heavy-ion collisions have recently been 
measured [HI EI], the need for extracting information that can be utilized 
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to probe the gravitational dual of the QGP has become pertinent. In the 
context of the gauge-gravity duality, the simplest field-theory objects are 
the connected correlation functions for the stress (energy-momentum) tensor. 
These are dual to graviton n-point functions j3]. 

We have recently calculated the one-particle irreducible (1PI) on-shell 
amplitudes of gravitons [TTJ and since extended the calculation to several 
connected n-point functions [12]. This work provides the basic ingredients to 
evaluate the connected stress-tensor correlators. We have focused on tensor- 
graviton amplitudes and on a certain kinematic region in which the am- 
plitudes simplify significantly. The calculation was performed for Einstein 
gravity and for its leading-order corrected theory, Gauss-Bonnet gravity. We 
have identified the following universal properties: 

1) All odd n-point functions of tensor gravitons vanish. 

2) For Einstein gravity, the four- and higher-point functions all go as the 
Mandelstam center-of-mass variable s. 

3) For the Gauss-Bonnet correction, the four-point function goes as s 2 and 
the higher-point functions go as s(s + v), where v is another independent 
angle, interpreted as a "generalized Mandelstam variable" . 

2 The gravitational dual 

We wish to explain briefly how on-shell graviton n-point functions are calcu- 
lated. The details can be found in [TTJ [12] ■ 

The gravitational Lagrangian is expanded out in numbers of gravitons, 
where a graviton represents a small perturbation of the metric from 
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its background value, — > g^ v + . The n order of the expansion 
can be used to read off the graviton 1PI n-point function. Since we are 
eventually interested in the dual stress-tensor correlation functions, the bulk 
n-point functions are evaluated at some large but finite radius. We calculate 
amputated amplitudes in momentum space so the n-point functions are just 
functions of momentum. 

In practice, one chooses the following ansatz for the gravitons: h^ v = 
4>(r) exp [iut — kz] . Here, z indicates the direction of graviton propagation 
along the brane (or along any other spacetime slice at constant radius) and r 
is the usual AdS radial coordinate, ranging from r = at the horizon of the 
black brane to infinity at the AdS boundary. The remaining two spacelike 
directions will be denoted by x and y. For future reference, a, b, ■ ■ ■ = 
{t, x, y, z} . We work in the radial gauge, for which h rr = h ar = for 
any choice of a. This allows the gravitons to be separated into three distinct 
classes [13]: tensors, vectors and scalars. In the following, we will only be 
interested in the tensors h xy . 

The results are presented below for Einstein gravity in units for which 
l/16nG^ is set to unity and s = —k^hiii is the standard Mandelstam 
variable. The external gravitons are symmetrized. 

We will be interested in the kinematic regime of "high momentum" in 
which the correlation functions of the stress-tensor simplify considerably. 
Also, the comparison between the leading Einstein result and possible cor- 
rections becomes simpler, as the number of derivatives in the interaction 
vertices is emphasized. In this regime, u, k and s are large enough so as to 
dominate over all radial derivatives, whether with respect to the perturbative 
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modes or the background. And so restricting to high momenta means that 
k 2 , u 2 , s ^> 1/L 2 , where L is the AdS curvature scale. 

The simplest example is the two-point function, lim (h xy (k)h xy (—k)) e = 



with all the odd-point functions vanishing. 

Rather than the above 1PI functions, we are really interested in con- 
nected functions of gravitons, as explained later. The connected functions 
can be calculated from the 1PI functions by evaluating tree-level Feynman 
diagrams, and these will have counterpart Witten diagrams on the boundary 
(see Fig. 1). However, we have found in [12] that, for Einstein gravity in the 
high-momentum regime, many of the connected functions are actually equal 
to the 1PI functions. We have checked that this persists up to the eight-point 
function and have strong indications that this phenomenon could persist for 
all the higher-point functions. 

In addition to explicit calculations showing the equality of the 1PI and 
connected amplitudes, the following general argument further supports their 
equality. The generating function of the connected amplitudes is the Legen- 
dre transform of the generating function of the 1PI amplitudes. In general, 
there is no reason for these functions coincide. Einstein gravity in the high 
momentum regime is, however, special. General covariance and the universal 
appearance of two derivatives combine to severely restrict the form of either 
generating function [TT| ITT], and so they can be expected to be in agreement. 

In any event, the difference between the connected and 1PI functions can 




(^) k 2 . And, in the high-momentum region, the higher-point functions for 
2n = 4, 6, 8 . . . are given by 




(1) 
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only be in their numerical coefficients, while the angular (in) dependence of 
the amplitudes remain linear in s. So that, 

/ T \ 4n_1 

lim((M 2n )con = a n - a, (2) 

with a 4 = 9/4 , a 6 = 75/8 and so on. 

3 From bulk to boundary 

We now want to translate the previous findings into statements about the 
dual strongly coupled fluid. 

Since the stress (energy-momentum) tensor T a b of the boundary field 
theory is dual to the bulk graviton h a b [3], the connected on-shell gravi- 
ton n-point functions can be used to calculate the (symmetrized) connected 
correlation function of the stress tensor, {T aibl T a2b2 . . .T anhn ) Co n ■ 

We need to apply the standard rules of holographic renormalization [T5J 
[161 [T7] to the bulk n-point functions. First, evaluate the bulk n-point func- 
tions at some large but finite radius r = tq . Second, multiply by an 
appropriate conformal factor Q q which is determined by the conformal di- 
mensions of the operators. Third, subtract off any divergences. One then 
only retains the part that survives in the r — > oo limit. 

The conformal factor for the stress-tensor correlation function can be 
deduced from the asymptotic form of the metric. For an AdS brane geometry, 
lim ds 2 = \r] a bdx a dx b ] + ^dr 2 , which tells us that Q — r/L . The 
power is q = A — 3 , where A = Aj is a sum over conformal dimensions 
and the subtraction of 3 comes from the dimension of the metric determinant. 
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The operators are the gravitons for which = 2 and the derivatives for 
which Ay = 1 • 

Following the described procedure, we find that all powers of r and L 
are stripped away from our previous expressions. What is left is a quantity 
that is finite and well defined at the boundary. Hence, the first and third 
steps turn out to be of no direct consequence. Subleading contributions to 
the metric could be important in principle but, in our case, only show up in 
terms that are asymptotically vanishing. 

The end result is that the field-theory correlation functions exhibit the 
same angular dependence as their bulk correspondents, 

((T X y) 2n )con — a 2nS , (3) 

with all odd-point functions vanishing. The coefficients equal to their 

bulk counterparts from Eq. ([2]). We emphasize that our results are only valid 
in the high-momentum region. The full 3-point function for the stress tensor 
does not vanish, as must follow from conformal symmetry. 

We remark that the 1PI functions for the stress-tensor are not very mean- 
ingful because of complications that arise in defining and separating the 1PI 
functions from the connected ones. (See [IB] for a recent explanation.) 

It is instructive to understand the kinematics of the high-momentum 
regime from the boundary point of view. The AdS scale L can be identi- 
fied with the spatial extent of the fluid [TjJJ [20], and so the meaning of the 
high-momentum region is that fc -1 , oj~ 1 and s -1 / 2 are smaller than the fluid 
size. Hence, one looks at correlation functions on scales less than the size 
of the fluid and does so by imposing an infrared cutoff on the momenta and 
frequencies. 
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To consistently apply the gauge-gravity duality in an AdS black hole (or 
brane) background requires that the horizon is larger than the AdS scale or, 
equivalently, hotter than the surrounding AdS space. In equation, nTL ^> 
1 . To understand the origin of this requirement, one considers geometries 
in which the boundary theory is defined on S 3 x R. In this case, L is also 
the radius of the boundary sphere, and so an extra dimensionful parameter 
is introduced which breaks conformal invariance spontaneously in the same 
way that the temperature does. In this context, the ratio ttTL corresponds 
to the ratio of the size of the boundary fluid to the thermal wavelength. 

Recall the high-momentum condition kL ^> 1 and, since we are also 
interested in the hydrodynamic limit, k/nT <C 1 . It follows that <C 
■^Y^ <C 1 and the dimensionless parameter nTL, which measures how much 
hotter the BH is than AdS space, also determines the range of momenta that 
occupy the high-momentum regime. 

4 Leading Correction 

In the context of the gauge-gravity duality, one neglects 1/N corrections, so 
that the leading correction to on-shell amplitudes comes from four-derivative 
corrections to the Einstein Lagrangian. The physical content of these correc- 
tions depends on the coefficient of the Riemann-tensor-squared term, as the 
other four- derivative terms (Ricci-tensor-squared and Ricci-scalar-squared) 
can be eliminated with a suitable choice of field redefinitions [2Tj |22| ITT] . 
Since our interest is in unitary theories whose equations of motion are at 
most second order in time derivatives, we use this freedom in field redefini- 
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Figure 1: Witten diagrams for stress-tensor correlation functions. The circle 
represents the AdS outer boundary, the dots depict interaction vertices, single 
lines correspond to undifferentiated gravitons and double lines correspond to 
differentiated gravitons. (a) A disconnected four-point function, (b) A 1PI 
four-point function at order e°. (c) A 1P1 six-point function at order e. (d) 
A 1PR connected six-point function at order e°. (e) A 1PR connected eight- 
point function at order e. (f) A 1PR connected ten-point function at order 
e 2 . 

tion to put the corrections in the Gauss-Bonnet form. 

Let e be a dimensionless parameter that measures the relative strength 
of the Gauss-Bonnet corrections compared to the leading Einstein result. 
Then e ~ lp/L 2 , where l v is the 5D Planck length. The parameter e can 
be determined by the coefficient of the Riemann-tensor-squared term in the 
Gauss-Bonnet combination. On the other hand, e also appears in the ratio 
of the shear viscosity to the entropy density - = [1 — 8e] [231 121] • On 
dimensional grounds, a term in the bulk Lagrangian with 2(m+l) derivatives 
should go as e m . A sensible perturbative expansion then requires that e < 1 . 
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Since Einstein and Gauss-Bonnet gravity both have a two-derivative lin- 
earized field equation, their real distinction comes from the non-linear inter- 
action terms. Einstein gravity has only two- derivative vertices, a property 
that greatly restricts and simplifies its higher-point IPI functions [T^, [TT]. In 
contrast, the four-point and e-order higher-point functions of Gauss-Bonnet 
gravity have four derivatives and, at higher orders in e, there can be many 
more. This distinction between different powers of e is illustrated in Fig. 1. 

We have found that the two extra derivatives in Gauss-Bonnet gravity 
lead to a more complicated angular dependence. However, the expansion 
to a given order in e is still universal. For instance, to first order in e, one 
obtains a specific angular dependence that is induced by the terms quartic 
in derivatives. Each additional order in e brings in added complexity. 

To first order in e, we have found that 

lim((M 2n )f// = lim((/^nf P7 

r— too t— >oo 

2 (2n\ r[n+§] f L\ 4n+1 . 

(4) 

where v = — k± ^2^=5 * s a "generalized" Mandelstam variable which 
accounts for the additional number of gravitons that are involved in the 
scattering process. For the four-point function, one should set v = . 

We have also found that, just like for Einstein gravity, all 1PR functions 
with a single internal line (e.g., Fig. le) vanish up to first order in e. Meaning 
that, at order e, the connected functions are equal to the IPI ones up to at 
least In = 6 . Then, because the e-order theory is similarly constrained by 
general covariance and a strict number of derivatives, this agreement could 
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very well persist for larger values of n. 

Going from the bulk to the boundary, we have to holographically renor- 
malize each of the two terms in Eq. separately because of their different 
r dependence. Then 

((T xy f n ) G c B on = {{T xy r) E Con + ec n s(s + v) , (5) 

such that v — when n = 2 and the odd-number functions are vanishing. 
For example, C4 = 3/4 , c% = 315/8 and so on. 

The relevance of the corrections to the interaction terms is determined by 
the value of L 2 k 2 e ~ k 2 l 2 , whereas the consistency of treating string theory 
effectively as a theory of gravity requires that k 2 l 2 < 1 . So that, if the mo- 
menta are "Planckian", then the corrections become substantial. But, as the 
hydrodynamic approximation requires that L 2 k 2 e <C (7iLT) 2 e , it is really 
the value of (tcLT) 2 € which determines whether the higher-point functions 
will be substantially corrected. If (7rLT) 2 e <C 1 , then the corrections will 
not appear in higher-point functions, only in the two-point function. This will 
be important when we discuss the possibility of experimentally determining 
the corrections. 

5 Contact with experiment 

The multi-point correlation functions of the stress (energy-momentum) ten- 
sor are the key to testing our ideas. Imagine a heavy-ion collision that creates 
a drop of QGP fluid. The drop expands and cools and then, after it freezes, 
a variety of decay products (mostly pions) stream outward. The energy, mo- 
mentum and the velocity of the decay products can be measured (more easily 
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for charged particles than for neutral ones), so that it is possible to measure 
the angular distribution of energy and momentum for a large number of 
particles emanating from the QGP fluid. The correlation functions of the 
stress-tensor components of the fluid well before freeze out can then be de- 
duced with the standard methods, in a way similar to evaluating multiplicity 
correlation functions (e.g., [25]). 

Let us now consider the relevant components of the stress tensor T xy = 
T yx ; these being dual to the tensor graviton modes. We are, assuming that 
the fluid flows in the z direction. Our objective is then to express these 
components in terms of simple observables. 

For a conformal boundary theory, the stress-tensor can be expressed to 
leading order as (e.g. 



lab - 3 



„pcpd / 2 

Au a u b + r] ab - 47r 2—-^ d c u d + d d u c - -rj cd rj cd d e u ( 



s txT \ 3 



(6) 



where u a is the four-velocity of the fluid, P ab = r] ab + u a u b projects vec- 
tors onto directions perpendicular to u a , p is the energy density, T is the 
temperature and rj/s is the ratio of shear viscosity to entropy density. To 
arrive at this form, we have used the relations p = p/3 for the pressure and 
s = (p+p)/T = (4/3)p/T . Indices are raised or lowered with the Minkowski 
metric r] ab . 

As the coordinates x and y are the spacelike directions perpendicular to 
u a , we can calculate T xy by disregarding the first two terms in and then 
eliminating the diagonal components; the latter by contracting the stress- 
tensor with the "binormal vector" e a b = {1 if a ^ b ; if a = b }. Then, 
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defining 47r2$ = T x y and a polarization tensor = ^e ab P° a P d b 



b ■ 



we have 



(7) 



where d a is a dimensionless spatial derivative, measured in units of the tem- 
perature d = -^fd . 

As the above equation makes clear, $ can be determined by measuring the 
angular distribution of the energies, momenta and velocities of the collision 
products. 

6 Scales, couplings and kinematics 

As we are interested in applying the gauge-gravity duality to calculate stress- 
tensor correlation functions, it is necessary to insure that our analysis con- 
forms with the duality's region of validity. 

The holographic dictionary tells us that the number of colors N is related 
to the 't hooft coupling A and the Yang-Mills coupling qym of the field theory 
as A = Qym-N ■ The first restriction comes from the large- N limit, for 
which N — > oo, qym such that A is finite and large. This means that 
1/N corrections are neglected while 1/A corrections are considered small but 
finite. But, for the QGP, A ~ 10 while N = 3 and there appears to be 
a conflict of interests. Nevertheless, as many features of the duality seem to 
be insensitive to this discrepancy, it is standard operating procedure to treat 
iV as much larger than A (see [27] for a recent discussion). We will adhere to 
this philosophy. 

A further constraint comes from the hydrodynamic approximation, which 
is an implicit consequence of treating the QGP as a fluid. In this case, 
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all momenta and frequencies are expected to be substantially smaller than 
the QGP temperature T (c/, the end of Section [3]). However, let us adopt 
the optimistic viewpoint that the hydrodynamic regime can be loosened to 
(k/nT) 2 , (lo/tiT) 2 < 1 and, within the context of our results, s/{nT) 2 < 1 . 
For the QGP at the LHC, T ~ 400 MeV and L ~ 7 fm , so that 
7rTL ~ 15 and the range in s is about 200. This means that the opportunity 
to distinguish different powers of s from data is open. 

As discussed previously, the value of the dimensionless parameter e(irLT) 2 
determines whether the corrections can be substantial in higher-point func- 
tions. For the QGP, using the above estimates and e = 1/8(1 — 4m]/ s) , 
we find that e(7rLT) 2 ~ 25(1 — Anrj/s) . Meaning that the higher-point 
functions are potentially sensitive to small corrections away from Einstein 
gravity. 

The disconnected correlation functions (ignored so far) can be used for 
additional consistency checks when comparing our predictions with experi- 
ment. Consider, for example, the disconnected part of the four-point function 
(Fig. la): (hhhh)D is = -(^(hh) 2 or ( ' hh ^ ma = 6(hh) . Clearly, one can 
use the value of the disconnected four-point function as an alternative means 
for ascertaining the two-point function and therefore rj/s. This method works 
also for higher-point disconnected functions. 

7 Conclusion 

We have calculated the n-point correlation functions for the stress tensor of 
a strongly coupled conformal fluid via the gauge-gravity duality. We have 
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explained schematically how to compare these results to statistical tests on 
the quark-gluon plasma at the LHC and what are the possible outcomes of 
such a comparison. 

Given that the experimental tests are indeed carried out, there are three 
distinct possibilities: (i) The gravitational dual is shown, within errors, to 
be in the universality class of Einstein gravity, (ii) the gravitational dual is 
shown to be in the universality class of a corrected Einstein gravity or (iii) 
the correlation functions of QGP do not follow from either. Each of the three 
possibilities would provide very interesting information about the nature of 
the QGP and about the existence of a gravity dual. The first outcome is, 
of course, our expectation. The second would be a surprise and a wonderful 
opportunity to study "Planckian" physics with the QGP. The third outcome 
would indicate that, if the the quark-gluon plasma has a gravitational dual, 
it has to be in a more complicated universality class, with the implication 
that the low value of rj/s is perhaps less significant in this context. 
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